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摘 要 


随 着 “工业 4.0” 和 “智能 制造 2025” 的 推广 ， 智 能 装备 正成 为 高 端 装备 未 
来 的 发 展 方向 。 在 这 个 过 程 中 ， 故 障 诊断 识别 系统 作为 一 项 重要 的 研究 领域 ， 可 
以 实时 抑制 机 械 的 不 平衡 振动 ,从 而 实现 装备 的 自主 健康 。 这 种 系统 不 仅 可 以 提 
高 设备 的 运行 效率 和 生产 能 力 ， 还 能 够 降低 设备 维修 的 成 本 和 风险 。 因 此 ， 研 究 
和 开发 这 种 系统 对 于 推动 智能 制造 和 高 端 装备 制造 的 发 展 具 有 重要 意义 。 

大 型 旋转 机 械 的 状态 监测 与 故障 诊断 是 现代 企业 生产 管理 的 必 备 手段 , 通过 
科学 的 监测 和 诊断 ， 可 以 实现 设备 的 高 效 、 安 全 、 可 靠 运行 ， 为 企业 的 可 持续 发 
展 提供 有 力 保障 。 本 课题 以 旋转 机 械 转 子 系统 典型 振动 故障 为 研究 主体 ,考虑 对 
其 进行 分 析 , 建立 旋转 机 械 转 子 系统 典型 故障 诊断 智能 识别 系统 。 本文 主要 研究 
内 容 与 结论 如 下 : 

《1) 本 文通 过 查阅 相关 文献 对 当前 的 机 械 智能 故障 诊断 算法 和 运用 神经 网 
络 进行 故障 分 类 的 研究 方法 进行 了 总 结 探 究 ， 并 最 终 决 定 选 择 RBF 神经 网 络 进 
行 故障 诊断 分 类 ; 

(2) 使 用 PCA 降 维 技术 对 提取 特征 后 的 数据 进行 降 维 处 理 , 来 降低 数据 的 
维度 来 解决 “ 维 数 灾难 ”问题 ， 使 数据 更 加 可 靠 ; 

(3) 通过 使 用 人 工 神经 网 络 中 的 径 向 基 神 经 网 络 建立 故障 分 类 软件 ， 利 用 
其 对 故障 通过 分 类 器 识别 并 进行 相关 实验 验证 ， 完 成 分 类 软件 设计 ; 

(4) 搭建 试验 台 采 集 数 据 进行 相关 实验 验证 ， 实 验 表 明 ， 诊 断 准确 率 达 到 
85% 以 上 。 


关键 词 : 典型 振动 故障 ; PCA 降 维 ; RBF 


第 1 章 绪论 


通过 状态 监测 与 故障 诊断 , 可 以 及 时 发 现 设 备 的 异常 运行 情况 ,避免 设备 因 
故障 造成 停机 维修 ， 从 而 减少 生产 损失 和 维修 成 本 以 及 提高 生产 效率 。 在 大 型 旋 
转机 械 的 状态 监测 与 故障 诊断 中 , 应 用 先进 的 传感器 、 信 号 处 理 和 数据 分 析 技 术 ， 
可 以 实现 对 设备 的 高 效 、 精 准 监 测 和 诊断 。 同 时 ， 还 可 以 通过 云 计 算 、 物 联网 等 
手段 ， 提 高 设备 的 运行 效率 。 总 之 ， 大 型 旋转 机 械 的 状态 监测 与 故障 诊断 是 现代 
企业 生产 管理 的 必 备 手段 , 通过 科学 的 监测 和 诊断 ,可 以 实现 设备 的 高 效 、 安 全 、 
可 靠 运行 ， 为 企业 的 可 持续 发 展 提供 有 力 保障 。 

见 代 工业 生产 自动 化 程度 高 ， 机 械 设 备 规模 、 结 构 和 工作 强度 逐渐 增 大 ， 长 


丙 
时 间 运 转 可 能 面临 复杂 多 样 的 故障 。 因 此 , 对 旋转 机 械 故 障 诊断 提出 更 高 的 要 求 。 

就 旋转 机 械 方面 的 故障 诊断 而 言 ， 国 际 上 已 经 有 了 相关 的 研究 成 果 。 例 如 ， 
美国 的 Sohre J. S. 发 表 论 文 总 结 了 旋转 机 械 的 典型 故障 特点 和 原因 ， 并 被 广泛 
应 用 于 高 速 旋转 机 械 故 障 诊断 领域 。 在 国内 方面 , 高 金吉 博士 深入 探索 了 高 速 旋 
转机 械 的 故障 机 制 , 并 基于 此 , 建立 起 一 种 以 10 类 58 种 为 基础 的 故障 类 型 的 分 
类 模型 ， 以 此 来 更 好 地 掌握 其 发 生 的 可 能 性 。 

关于 旋转 设备 故障 诊断 的 方法 , 主要 包括 如 下 几 种 : 常规 诊断 法 、 信 息 
方法 、 采 用 专家 系统 、 采 用 神经 网 络 、 采 用 模糊 方法 和 可 能 性 理论 的 方法 等 。 这 
些 方法 可 以 帮助 我 们 更 准确 地 诊断 机 械 故 障 ， 提 高 故障 处 理 效 率 和 准确 性 。 

转子 不 平衡 是 旋转 机 械 振动 的 一 个 重要 影响 因素 , 特别 是 对 于 高 速 旋转 机 械 。 
同时 ,转子 不 平衡 也 会 影响 其 它 转子 故障 的 发 生 和 加 剧 。 为 了 解决 这 一 问题 , 现 
有 在 线 自动 平衡 系统 可 以 实时 抑制 机 床 不 平衡 振动 ， 采 用 电磁 式 自动 平衡 方式 ， 
能 够 在 短 时 间 内 抑制 机 床 主轴 的 不 平衡 振动 。 这 种 自动 平衡 系统 可 以 通过 测量 机 
床 的 振动 信号 ， 计 算出 机 床 的 不 平衡 情况 ,并 通过 控制 电磁 力 的 大 小 和 方向 来 实 
现 对 机 床 的 自动 平衡 。 通 过 引入 先进 的 自动 平衡 技术 , 我 们 能 够 大 幅 改 善 机 器 的 
运行 状态 , 并 且 能 够 更 准确 地 控制 零件 的 位 置 。 这 种 技术 不 仅 能 够 更 好 地 保证 设 
备 的 正常 运行 ， 还 能 够 大 幅 提 升 生产 效 率 ， 带 来 更 多 的 经 济 收 入 。 

本 课题 拟 对 旋转 机 械 转 子 系统 典型 故障 振动 系统 进行 开发 , 通过 研究 降 维 方 
法 的 学 习 使 用 ， 神 经 网 络 的 构建 与 验证 提出 旋转 机 械 转 子 系统 故障 诊断 的 方案 ， 
以 便 提高 旋转 机 械 的 工作 效率 ， 保 障 其 安全 可 靠 运行 。 


第 2 章 RBF 和 了 PCA 理论 基础 


2.1 PCA 降 维 

首先 介绍 下 为 什么 要 进行 数据 降 维 : 

通过 对 高 维 数据 的 分 解 ， 可 以 使 其 变 得 更 加 简单 、 清 晰 。 这 种 方法 的 优点 包 
括 : (1) 节约 运算 成 本 费用 : 分 解 高 维 数据 所 消耗 的 存储 空间 和 运算 资源 较 小 ， 
从 而 节约 了 运算 成 本 费用 ; (2)〉 改善 模型 性 能 : 分 解 高 维 数据 使 得 模型 变 得 简 
单 ， 从 而 避免 了 过 拟 合 的 问题 ，(3) 去 除 见 余 信 息 : 通过 采取 数据 降 维 技术 ， 
我 们 能 够 消除 多 维度 的 数据 ,从 而 改善 模型 的 精度 与 可 靠 性 ，(4) 可视化 展示 : 
高 维 数据 难以 进行 可 视 化 展示 ， 而 低 维 数据 可 以 更 容易 地 进行 可 视 化 展示 ， 有 助 
于 人 们 更 好 地 理解 和 分 析 数 据 ; (5) 数据 压缩 : 通过 采取 数据 压缩 技术 ， 将 多 
维度 的 数据 维度 降低 ， 从 而 节约 了 存储 空间 。 

PCA(Principal Component Analysis, 主 成 分 分 析 ) 作 为 一 种 高 级 的 降 维 算法 ， 
它 可 以 有 效 的 把 复杂 的 高 维 数据 分 解 成 更 加 精细 的 维度 , 从 而 更 加 精准 的 抓 取 出 
数据 的 关键 特征 。 PCA 算法 旨 在 通过 线性 变换 , 把 原 有 的 数据 投影 至 一 个 全 局 的 
位 置 ， 以 达到 减少 其 变形 量 的 目的 ,并且 可 以 有 效 地 提高 投影 的 精度 ， 以 达到 有 
效 的 数据 降 维 。 PCA 算法 的 优点 是 可 以 消除 数据 中 的 元 余 信 息 和 噪声 , 提高 模型 
的 准确 性 和 解释 性 。 同 时 ，PCA 算法 也 可 以 用 于 数据 压缩 和 可 视 化 等 方面 的 应 
用 。 

PCA 降 维 的 基本 原理 如 下 图 2-1 所 示 : 

首先 对 原始 数据 进行 去 均值 化 操作 ,使 得 数据 中 心 为 原点 。 经 过 计算 ， 首先 
可 以 得 出 协 方差 矩阵 ， 它 描述 了 各 个 维度 之 间 的 相互 关系 。 通 过 将 协 方差 窍 阵 的 
特征 值 进行 分 解 ,就 可 以 更 加 清晰 地 了 解 每 一 个 特征 值 和 它们 的 相关 向 量 。 为 了 
更 好 地 表示 原始 数据 , 可 以 考虑 将 其 投影 到 一 个 新 的 矩阵 上 , 从 而 实现 降 维 处 理 。 
在 PCA 算法 中 , 特征 向 量 是 PCA 算法 的 核心 , 它 代表 了 数据 在 新 坐标 系 下 的 方 
向 。 特 征 向 量 的 长 度 代 表 了 数据 的 方差 ,方向 则 代表 了 数据 的 主要 变化 方向 。 通 
过 PCA 算 法， 就 可 以 实现 从 原始 数据 中 选择 前 K 个 最 大 的 特征 值 ， 并 将 其 转换 成 
具有 较 高 方差 的 特征 向 量 ， 从 而 有 效 地 实现 数据 的 降 维 。 这 样 ， 最 后 就 可 以 将 具 
有 较 高 方差 的 特征 维度 保留 下 来 , 而 将 具有 较 低 方 差 的 特征 维度 排除 在 外 ， 从 而 
达到 降 维 的 目的 。 


图 2-1 PCA 降 维 原 理 图 


实现 PCA 的 一 种 算法 如 下 : 
输入 :数据 集 X ={x,x%,%,.…,x%}， 需 要 降 到 k 维 。 
(1) 去 平均 值 ， 即 所 有 特征 减 去 元 。 
(2) 计算 协 方差 矩阵 二 XX” , 这 里 处 理 样本 数量 n 或 n 一 1 与 否 ,与 特征 向 
量 无 关 。 
(3) 解决 协 方差 矩阵 - XX7 的 关键 问题 ， 包括 确定 其 特征 值 和 方向 。 


(4) 特征 值 排序 , 将 对 应 的 按 大 小 顺序 排列 的 前 k 个 特征 向 量 按 行 排列 

为 特征 向 量 和 矩阵 P。 

(5) 通过 将 k 条 特征 值 映射 至 一 组 新 的 空间 ， 使 得 Y=PX 成 为 可 能 。 
实现 步骤 画 出 流程 图 如 下 : 


减 去 平均 值 (中心 化 ) 


计算 协 方差 矩阵 


v 


计算 特征 值 和 特征 向 量 


将 特征 值 排序 


保留 前 N 个 最 大 的 特征 值 对 应 的 特征 向 量 


将 原始 特征 转换 到 N 个 特征 向 量 构建 的 新 空间 


图 2-2 PCA 降 维 实现 算法 


其 中 ， 样 本 均值 : 


样本 方差 : 
Be 
S i x) 
样本 X 和 样本 Y 的 协 方差 
cov(X,Y) =E[(X—E(X))(Y -EY))] 
1] 之 a i 
2 


2.2 径 向 基 神 经 网 络 (RBF) 

RBF 神经 网 络 的 架构 包含 三 个 基本 单元 : 输入 、 处 理 器 、 以 及 反馈 。 在 处 理 
器 的 基本 单元 中 ,输入 是 指 从 系统 获取 的 数据 ， 它 会 被 转换 成 一 系列 的 模型 ， 然 
后 被 处 理 器 处 理 。 在 处 理 器 的 基本 单元 (Radial Basis Function〉 中 ， 它 们 会 使 用 
径 癌 基 函 数 来 处 理 高 维 数 据 。 径 同 基 函数 是 一 类 常见 的 数学 方法 , 其 中 包含 高 斯 
函数 和 多 项 式 函 数 。 反 馈 层 的 节点 数 通常 与 目标 变量 的 维度 相同 , 用 于 输出 网 络 
的 预测 结果 。 反馈 层 的 节点 通常 采用 线性 函数 作为 激活 函数 , 将 处 理 层 的 输出 进 
行 线性 组 合 , 得 到 最 终 的 预测 结果 。RBF 神经 网 络 的 显著 特性 在 于 它 的 模型 准确 
性 和 高 效 性 , 它 的 模型 能 够 很 好 地 模拟 任何 连续 函数 , 并 且 在 短 时 间 内 就 能 获得 
准确 的 结果 ， 避免 了 模型 的 局 限 性 。 另 外 ， 它 还 适合 在 复杂 的 情况 下 进行 模型 的 
训练 ;， 不仅 如 此 ， 对 于 多 分 类 问题 ， 可 以 使 用 一 次 训练 多 个 输出 。 

RBF 神经 网 络 结构 如 图 2-3 所 示 。 
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图 2-3 RBF 神经 网 络 结构 


当 输 入 向 量 瑟 =[zz zx], 输出 向 量 了 =[y,,y,,…,y] 时 , RBF 神经 网 络 


的 隐 含 屋 和 结果 为 : 


yj) => oD(x) jel,2,..,l 


k=1 


D, (2X) 一 en- 


20 


其 中 : o 为 基 函 数 宽度 cy 为 隐 含 层 第 有 个 神经 元 与 输出 层 第 7 个 神经 元 之 
间 的 连接 权 值 ; c, 为 第 大 个 神经 元 的 中 心 点 向 量 ，@, (0 为 第 严 个 径 向 基 函 数 的 
Gauss 函数 。 


第 3 章 基于 PCA-RBF 的 代码 实现 与 实验 验证 


3.1 基于 MATLAB 的 PCA 特征 降 维 
在 MATLAB 中 可 以 通过 编写 上 述 算法 代码 来 实现 PCA 数据 降 维 ， 如 下 图 
3-1 所 示 : 


function S=princa(X) 
[m,n]=size(X); gp ey 

%------------- : 标准 化 矩阵----------------- % 
mv=mean(X); x 计算 各 变量 的 均值 
st=std(X); % 计 算 各 变量 的 标准 差 
X=X-repmat(mv,m,1); % 标 准 化 矩阵 X 

全 第 二 步 : 计算 相关 系数 矩阵 ----------------- % 
% R1=X'*X/(m-1); % 方 法 一 : 协 方差 矩阵 计算 公式 
% R2=cov(X); % 方 法 二 : 协 方差 矩 阵 计算 函数 
R=corrcoef(X); % 方 法 三 : 相关 系数 矩阵 函数 
%------------- 第 三 步 : 计算 特征 向 量 和 特征 值 ----------------- 
[V,D]=eig(R); % 计 算 拭 阵 R 的 特征 向 量 粉 阵 V 和 特征 值 矩 阵 D， 特征 值 由 小 到 大 

% 将 特征 向 量 拭 阵 V 从 太 到 小 排序 
% 将 特征 值 矩 阵 由 大 到 小 排序 


E=diag(D); % 将 特征 值 矩 阵 转 换 为 特征 值 向 量 
ee 第 四 步 : 计算 贡献 率 和 累计 贡献 率 ----------------- % 
ratio=8; % 累 计 贡 献 率 

for k=1:n 


r=E(k)/sum(E);  ”% 第 k 主 成 份 贡 献 率 
ratio=ratio+r; %% 索 计 贡 献 率 
if(ratio>=8.8) % 取 累计 贡献 率 太 于 等 于 96% 的 主 成 分 
q=k; 
break; 


i 第 五 步 ; 计算 得 芬 ===== ==- % 
V=V(:,1:q); 
S=X*V; 

图 3-1 PCA 降 维 算法 MATLAB 实现 


在 MATLAB 中 也 可 直接 调用 pca 函数 实现 PCA 数据 降 维 ， 其 调用 格式 如 


[coef!, SCOre, latent| = pca (X) 


降 维 效果 无 明显 差异 ， 本 课题 选用 调用 pca 函数 实现 PCA 数据 降 维 。 
3.2 基于 MATLAB 的 径 向 基 神 经 网 络 故障 分 类 实例 


在 MATLAB 中 男 出 RBF 径 向 基 传 递 函 数 如 下 图 3-2 所 示 : 


径 向 基 传 递 函数 


传递 函数 
图 3-2 径 向 基 传 递 函数 〈Gauss 函数 ) 
RBF 神经 网 络 的 训练 过 程 主 要 包括 以 下 几 个 步骤 : 中 确定 隐 含 层 节 点 数 : 通 
过 交叉 验证 等 方法 ， 确 定 隐 含 层 节 点 数 ， 以 提高 神经 网 络 的 性 能 和 泛 化 能 力 。 包 
初始 化 权重 : 随机 初始 化 隐 含 层 和 输出 层 的 权重 ,以 便 开始 训练 神经 网 络 。@@® 计 
算 径 向 基 函 数 : 通过 非 线 性 变换 处 理 输入 数据 ， 计 算 径 向 基 函 数 的 值 ， 以 获得 更 
准确 的 结果 。 外 计算 输出 层 权 重 : 通过 最 小 二 乘法 等 方法 ， 计 算 输 出 层 的 权重 ， 
以 产生 最 终 的 输出 结果 。@@ 更 新 权重 : 通过 反 向 传播 算法 等 方法 ， 更 新 神经 网 络 
中 的 权重 ， 以 提高 神经 网 络 的 性 能 和 准确 性 。 
在 MATLAB 中 可 调用 函数 newrp 实现 RBF， 其 调用 格式 如 下 : 
net =newrb(P,T, goal, spread, MN, DF) 


es 


式 中 部 分 参数 解释 如 下 表 所 示 : 


表 3-1 RBF 参数 解释 


参数 实际 意义 

P 输入 向 量 

克 目标 分 类 向 量 

goal 均 方 误差 目标 ， 默 认为 0 
spread 径 向 基 函 数 的 扩展 速度 ， 默 认为 1 
MN 径 向 基 神 经 元 最 大 值 

Dr 径 向 基 增 加 步 长 ， 默 认为 25 
net 返回 值 ， 即 一 个 RBF 网 络 


通过 构建 3 维 数据 集 , 我 们 从 300 组 不 同 的 不 平衡 和 不 对 中 的 情况 下 收集 各 
300 组 信息 ,以 及 300 组 的 干扰 , 构建 出 P=900x3 维 的 模型 ， 以 期 望 能 够 更 好 地 
模拟 故障 诊断 情况 。 将 P 输入 集 随机 分 成 540 组 和 360 组 ， 分 别 用 于 进行 训练 
和 评估 。 

RBF 神经 网 络 分 类 部 分 使 用 MATLAB 编写 的 部 分 代码 如 下 : 


% 2. 随机 产生 训练 集 和 测试 集 

temp = randperm(size(feature,1)); 

% 训练 集 一 999*6.6 个 样本 

P_train = feature(temp(1:540),:)"; %P 和 T 列 代表 个 数 必 须 保持 相等 
T train = labal(temp(1:540),:)"; 

% 测试 集 -366 个 样本 

P_test = feature(temp(541:end),:)'; 

T_test = labal(temp(541:end),:)'; 

N = size(P_test,2); 


‰yV III.RBF 神 经 网 络 创 建 及 仿真 测试 

X% 

% 1. 创建 网 络 

t = cputime; 

rbf_net=newrb(P_train,T train,8.864,0.82,158,1); 


X% 

% 2. 仿真 测试 

T_sim = sim(rbf_net,P_test); 

result=round(T_sim); 

t = cputime - t; 

党 TV 性 能 评 众 

X% 

% 1. 相对 误差 error 

% 计 数 判 断 正 确 则 为 1， 否 则 为 8 

no_test=size(result,2); 

count=0; 

for n=1:1:no_test 

if result(:,n)==T_ test(:,n) 
count=count+1; 

else 


end 
end 
accuracy=count/no test; 


图 3-3 RBF 神经 网 络 分 类 MATLAB 实现 
在 MATLAB 中 对 其 进行 训练 ， 如 下 图 3-4 所 示 : 
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3-4 RBFNN 的 训练 过 程 
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3.3 实验 验证 


把 采集 的 不 平衡 状态 、 不 对 中 状态 下 数据 结合 之 前 编 的 算法 进行 算法 验证 。 

采用 上 述 典 型 的 振动 智能 诊断 算法 , 可 以 有 效 地 检测 出 不 平衡 和 不 对 中 的 数 
据 ,， 划 分 训练 集 测 试 集 之 后 ,从 中 抽取 各 80 组 数据 , 不 对 中 故障 标签 设置 为 1， 
不 平衡 故障 标签 设置 为 2， 并 将 其 结果 展示 在 图 3-5 中 ， 从 而 更 加 清晰 地 了 解 它 
们 的 状态 ， 如 下 图 所 示 : 


图 3-5 不 对 中 (上 )〉 和 不 平衡 (下 〉 故障 诊断 识别 结果 
经 过 多 次 实验 验证 分 析 , 可 以 发 现 典型 的 振动 智能 诊断 软件 能 够 准确 地 识别 
出 不 对 中 和 不 平衡 的 故障 情况 ,其 平均 识别 率 达 到 85% 以 上 , 完全 符合 课题 的 预 
期 要 求 。 


第 4 章 结论 


本 文 以 旋转 机 械 转子 系统 典型 振动 故障 为 研究 主体 ， 从 理论 、 仿真 和 试验 等 
多 个 方面 建立 旋转 机 械 转 子 系统 典型 故障 诊断 智能 识别 系统 并 进行 相关 实验 验 
证 ， 并 取得 了 良好 的 实验 效果 。 

对 提取 出 的 特征 进行 PCA 降 维 处 理 等 操作 ， 并 在 MATLAB 软件 中 编写 相 
关 代码 实现 这 一 功能 ; 在 完成 特征 降 维 之 后 , 需要 选择 合适 的 分 类 器 进行 故障 分 
类 。 本 文通 过 对 人 工 神经 网 络 中 的 多 种 分 类 器 网 络 进行 网 络 结构 探究 、 隐 含 层 传 
递 函数 比较 等 进行 多 种 人 工 神 经 网 络 的 优 缺 点 对 比分 析 , 并 最 终 选 择 径 向 基 神 经 
网 络 (RBF) 作为 该 课题 的 故障 分 类 器 。 之 后 对 降 维 后 的 数据 进行 神经 网 络 的 训 
练 与 验证 ， 并 在 MATLAB 中 编写 相关 代码 完成 这 一 功能 。 

最 后 ， 为 验证 算法 的 效率 和 准确 性 以 及 集成 方案 的 可 行 性 ,考虑 做 相关 实验 
进行 实验 验证 。 将 不 平衡 、 不 对 中 这 两 种 情况 下 的 故障 数据 进行 采集 ， 然 后 ， 带 
入 到 原来 的 故障 诊断 算法 中 进行 验证 , 故障 诊断 率 达 到 85% 以 上 ,符合 预期 要 求 。 


